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Abstract 

We derive the axial Ward identity for lattice QCD with domain- wall fermions, and from which we 
obtain a formula for the residual mass (|45p -(|46 p . that can be used to measure the chiral symmetry 
breaking due to the finite extension N s in the fifth dimension. Furthermore, we obtain an upper 
bound for the residual mass in lattice QCD with the optimal domain-wall fermion. 
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I. INTRODUCTION 



The chiral symmetry of massless fermion field plays an important role in particle physics. 
It forbids the additive mass renormalization which causes the fine-tuning problem associated 
with the scalar field. In QCD, the chiral symmetry [SU^Nf) x SUn(Nf)] of Nf massless 
quarks is spontaneously broken to SUv{Nf), due to the strong interaction between quarks 
and gluons. This gives the (nearly) massless Goldstone bosons (pions) and their specific in- 
teractions. To investigate the spontaneously chiral symmetry breaking (or hadronic physics) 
in QCD, it requires nonperturbative methods. So far, lattice QCD is the most promising 
approach. However, in lattice QCD, formulating lattice fermion with exact chiral symmetry 



at finite lattice spacing is rather nontrivia 



(DWF) on the (4+l)-dimensional lattice and the overlap fermion on the 4-dimensional 
lattice Q. 

For lattice QCD with DWF, in practice, one can only use a finite number N s of sites 
in the fifth dimension. Thus the chiral symmetry of the massless quark fields is broken, 
and the emergent question is whether the chiral symmetry is preserved optimally. The 
answer is negative since the effective 4-dimensional Dirac operator of the conventional DWF 
corresponds to the overlap Dirac operator with the polar approximation of the sign function 
of if. 

In 2002, one of us (TWC) constructed the optimal domain-wall fermion (ODWF) {^J 
such that the effective 4D lattice Dirac operator attains the mathematically optimal chiral 
symmetry for any finite N s , exponentially- local for sufficiently smooth gauge backgrounds 
4j, and independent of the lattice spacing in the fifth dimension. The basic idea of ODWF 
is to construct a set of analytical weights, {u s , s = 1, ■ • • , N s }, one for each layer in the fifth 
dimension, such that the chiral symmetry breaking due to finite N s can be reduced to the 
minimum. The 4-dimensional effective Dirac operator of massless ODWF is 

D = ^-[l+ l5 S opt (H)}, 

oA J ~i+nr=iT,' s "i+^' 

which is exactly equal to the Zolotarev optimal rational approximation of the overlap Dirac 
operator. That is, S opt (H) = HR Z (H), where Rz(H) is the optimal rational approximation 

of (H 2 y i/2 B,B- 



which is realized by the domain-wall fermion 
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However, in the original ODWF formulation (3j, the valence quark propagator cannot 
be expressed in terms of the correlation function of the quark fields defined in terms of the 
boundary modes, unlike the conventional domain-wall fermion. In 2003, one of us (TWC) 
solved this problem by introduced two transparent layers with u s = |7j, as boundary layers 
appending to the original action of ODWF such that the quark fields defined in terms of these 
two transparent layers obey the usual chiral projection rule in the continuum, independent 
of the gauge fields. Consequently any observable constructed with the quark fields manifests 
the symmetries exactly as those of its counterpart in the continuum. The salient feature 
of a transparent layer (with u s = 0, and T s — 1) is that its presence does not change the 
effective 4D Dirac operator. 

In this paper, we derive the axial Ward identity for lattice QCD with ODWF. We find 



that it is necessary to extend the idea of transparent layers introduced in Ref. to add 
another two transparent layers at the central region of the fifth dimension. With these four 
transparent layers, the action of lattice QCD with ODWF can be written as 

N s +3 

= S ^4>s{x){{p s D w + T) X , X 'S S , S > + (cr s D w - l) XtX '(P-5 8 f jS+1 + P+5 s / )S _i)}'0 s '(^ / )) W 

s,s'=Q x,x' 

with boundary conditions 

P+ip(x, -1) = ~rm q P + ip(x, N s + 3), 

(2) 

P-ip(x, N s + 4) = -rm q P^ip(x, 0), 

where P± = (1 ±7s)/2, D w is the standard Wilson-Dirac operator plus a negative parameter 
— m (0 < m < 2), m q is the bare quark mass, and r is a parameter depending on {p s , a s } 
and m such that the valence quark propagator agrees with (7 M d M + m q )~ l in the continuum 
limit. The two transparent layers at the boundaries are specified by imposing p = Pn s +3 = 
cr = 0"tv s +3 = 0. The two additional transparent layers can be located at s = n, and 
s = n + 1, where n = [N s /2]. In other words, they have p n = p n+ \ — a n — cr„ +1 = 0. In the 
original ODWF formulation [3J, the nonzero p s and a s are set to be u s (the optimal weight). 
The quark fields are defined in terms of the boundary modes 

q(x) = [P-ipo{x) + P+^n s +3(x)} , 
q(x) = y/r [ip (x)P+ + ipN s +3(x)P_] . 
Following the derivation given in Ref. , it is straightforward to show (in Section IHIj) 



that the valence quark propagator in a gauge background is equal to the correlation function 



of the quark fields, i.e., 

(q(x)q(y)) = (D c + m q )-\x, y), (4) 



where 



r 1 - y 5 S 



1 TtV s +3 rp 

S = lls =° > (6) 
1 — H. 



H s = (p s + (7 s )H w [2 + (p s -a s )j 5 H w ] \ H w = y 5 D w . (8) 

Obviously, a transparent layer (with p s = u s = 0) does not change 5 and D c since its T s = f . 
Setting the nonzero weights p s = c u s + d, and a s = c u s — d, where c and d are constants, 
then if, = u s H, 

H = cH w (l + d l5 H w )~\ (9) 

and the parameter r entering the boundary conditions (j2J) is fixed to r = [2?71q(1 — dmo)]^ 1 
such that (_D C + in the free fermion limit agrees with (7 M <9 M + rrig)^ 1 in the continuum 

limit. Moreover, for H s = u s H, interchanging any two layers in the fifth dimension gives 
the same S, since {T s j commute among themselves. 

n 

For finite N s , with the optimal weights {u s } given in Ref. [3], S is exactly equal to 
the Zolotarev optimal rational approximation of the sign function of H, i.a, S = S op t(H) = 
HRz(H), where Rz(H) is the optimal rational approximation of (if 2 ) -1 / 2 5|,|6]- In the limit 
N s — > oo, S — > H(H 2 )~ l l 2 , and D c becomes exactly chirally symmetric, and (D c + m q )~ l is 
well-defined for nonzero m q , even though D c is ill-defined for topologically nontrivial gauge 
background [8j. 

In practice, only the case d = gives if = cH w (without the denominator), which is 
much easier for the projection of the low-lying eigenmodes of D — D c (l + rD c )~ 1 than other 
cases with d ^ 0. Since the low-lying eigenmodes of D are vital for extracting many physical 
observables, the original formulation 3| with d = (and c = 1) seems to be a good choice. 

We note in passing that setting the nonzero weights p s = c\ (constant) and <j s = 
C2 (constant) covers all variants of conventional do main- wall fermions, with S equal to the 
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polar approximation of the sign function of H, 



where 



S polar (H) 



tt 2 V^n 



H 2 +d[ ' 



N s = 2n + 1 (odd) 
N* = 2n (even) , 



bi = sec 2 



7T 



dj = tan z 



7T 



l- 



For example, setting p s = 1 and cr s = 0, ([T]) reduces to the conventional domain- wall fermion 
with H = H w (2 + 75ifu,) _1 jjj, and p s = <r s = 1 to the Borici's variant with H = H in flO ]. 
and p s = c + d and <r s = c — d to the Mobius variant with H = cH w (l + d / -f^,H w )~ 1 11]. 



II. AXIAL WARD IDENTITY 

Now we consider Nf flavors of quarks with degenerate mass m q , and the infinitesimal 
flavor non-singlet transformation 

5ip s (x) = i9 s (x)\ a ip s (x), 

(10) 

8ip s (x) = -ie s {x)ip s {x)\ a , 



where 



e s {x) 



6(x), < s < n = , 
-9(x), n + l<s<N s + 3. 



Here A a is one of the flavor group generators in the fundamental representation, and the 
flavor indices of ip s (x) and ips{x) are suppressed. Under the transformation fflOl) . the change 
of the action consists of the following three parts: 

N s +3 

s E X> [ux)x a D w (x, y )Uy)] = 

s=0 x,y x fi 

N s +3 



5 l^s(x)\ a P-^ s+1 (x) - ^ s (x)A a P + V> s -i(x)] = - J22i9(xM(x,n) + m q q(x)X a l5 q(x)}, 

s=0 x 
N s +3 

5 E E^^( x ) A " jD -( a; '^[ p -^+ i ^) +p +^- i ^)]> = E^( x )E a mW' ( n ' 



s=0 
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where 



Aa»/(:e, s) = /(x, s) - /(a; - /i, s) 

N s +2 



sign - s + - J s] 



J 5 a (x, n) = -Mx)\ a P^ n+1 (x) + 4> n+1 {x)\ a P + ^ n {x), (12) 

k;(x) = k; + (x) + k a -(x), 



Ns+2 , 



kf(x, s) = y [^(^)A a (l - 7m)^(^)^ sT i(x + //) - + /i)A a (l + 7 , t )^(x)P±^ Tl (x)] . 

Now the role of the two transparent layers at s = n and s = n + 1 becomes obvious. 
If we want to keep Jg (I12p not depending on D w (similar to the Jg in the conventional 
DWF) and to express (II II) in terms of the divergence of a 4-current, then it is inevitable 
to introduce two transparent layers in the central region of the 5th dimension. This can be 
seen as follows. For 1 < s < n — 1, or n + 2 < s < iV s + 2, we have 

5^2a s tp s (x)X a D w (x,y)P-ilj s+1 (y) = Ti^-4) 8 {x)\ a a s [0{x)D w {x,y) - D w (x, y)9(y)]P-ip a+1 (y), 

x,y x,y 

which can be written in the form of Y2 X 9{x)A^J^. However, at s — n, it gives 
5^2cr n 'ip n (x)\ a D w (x,y)P-'i/j n+1 (y) = -i^2^ s (x)\ a a n [9(x)D w (x,y) + D w (x,y)6(y)}P-ip n+1 (y), 

x,y x,y 

which cannot be expressed in terms of the divergence of a 4-current unless o n = 0. Similarly, 
we also set a n+ \ = 0. Furthermore, for consistency, we must also set p n = p n+1 = such 
that T n = T n+1 = 1. 

For any observable O, the variation of its vacuum expectation value with respect to (flOl) 
must vanish, i.e., 5 a (0) = 0, which gives the axial Ward identity 

A„< JZ(x)0(y)) = 2mM^ a l^)0(y)} + 2(J^(x, n)0{y)) + i(5 a O(y)}. (13) 

where J%(x) = k"(x) + j^(x). As N s — > oo, the anomalous term (Jg(x,n)0(y)) vanishes if 
0(y) only involves the quark fields, following the same argument given in Ref. [l^ . 
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After summing over all sites x, the LHS of (113]) vanishes, and its RHS gives 

-i J2(5 a O(y)) = 2m q ^(g(*)A a 7 5 g(*)(%)> + 2 ^(J 5 a (x, n)0{y)). (14) 

Thus, the effect of chiral symmetry breaking due to finite N s can be regarded as an additive 
mass to the bare quark mass m q , the so-called residual mass 

mM] = (15) 

which serves as a measure of the chiral symmetry breaking due to finite N s . In the limit 
N s — > oo, S(H) = H/ a/ H 2 and m res = 0. Obviously, in a gauge background, the residual 
mass ( TT5l) depends on the observable O as well as its location y. Thus it is necessary to take 
into account of the residual mass at all locations. This can be accomplished by summing 
over all lattice sites y in the axial Ward identity (fT^j) to obtain the global residual mass 

M , n] Y, x , y (Ji(x,n)Q(y)) 
^ 1 E,>>)>75^)£%)>- 1 j 

For 0(y) = q(y)X b j5q(y), ( Il5l) and ( fT6l) become 

/.a E x (J5(x,n)q{y)X b 'y 5 q(y)) 

[V) E, (q(x) X^ 5 q(x)q(y) X b l5 q(y)y 1 <} 

which are usually used as a measure of the chiral symmetry breaking due to finite N s . In the 
following, we will restrict our discussions to the residual mass ( IT?]) , and the global residual 
mass flUD. 

III. GENERATING FUNCTIONAL FOR n-POINT GREEN'S FUNCTION 

In order to express the residual mass ( TTTf) in terms of the quark propagator, we first derive 
the generating functional for the n-point Green's function of the fermion fields for lattice 
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QCD with ODWF. With four transparent layers, the action ([I]) can be rewritten as 

N s +3 

A f = $s( x h5{(p s H w P + + p s H w P- + P+- P-) x , x ,^ s (x) 

s=0 x,x' 

+(a s H w P„ + a s H w P + + P_ - P + ) x , x ,(P^ s+1 (x') + P+^tix'))} 
= ^^(^)7 5 {(Ps^P- + <? S H W P + - l) x ^[P + ^ s ^{x') + P^ s {x')\ 

s=0 x,x' 

+ (PsH w P + + a s H w P- + l) x , x >[P+Mx') + P-i> s+ i(x')}} 

s=0 x,x' 

where 

Q s ± = p s H w P± + a s H w P T ± 1. (19) 

Defining 

7] s = {P-5 S > !S + P + 5 s > >8 -i)ip 8 > ips = {P-S s ',s + P+S s ', s+ i)r]s', (20) 
Vs = tpslsQ- <^^ 8 = VsiQ-y 1 !^ (21) 
T s = -(Q;)- 1 ^ = H s = (p s + a s )H w [2 + (p, - a s ) lh H w }-\ (22) 

then the action ([1]) can be expressed in terms of rj, fj fields 

N 3 +2 

Af = ffo(P- - rm q P + )r] - f] r) + ^ IVsVs ~ fj s T~\ +l ] 

s=l 

+Vn s +sVn s +3 - rj Ns+3 (P + - rm q P^)r] , (23) 

where the space-time indices have been suppressed. 

In order to evaluate the Green's function of the fermion fields in the expression of the 
residual mass, we need to add the following external source terms to (|2"31 

VnJ-n + Jn+lVn+l + Vn+lJn+1, 

J q Q + qJ q = JVo ~ VoP+J + Vn s +3P-J, 

where 

J = y/rJ q , 

J = y/rJ q . 
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Then the generating functional for n-point Green's function is defined as 



Z[J q , J g , J n , J n+ i, J n+ i] = J I [drf\[drf\e Sj 



where 



Sj = Aj - Jr/o + T] P + J - r] n J n - J n+1 r] n+1 - r] n+1 J n+1 - r/ Ns+3 P^J, 
and J is the Jacobian of the transformation, 

Af s +3 

J= J] det(p s ^P_ + a s H w P + - 1). 

s=0 

Now using the Grassman integration formula 



J dxd X e-* Mx+ * x+ * v = e" M ' lv det M, 



and integrating (t] s ,T] s ) successively from s = N s + 3 to s — 1, (1241 becomes 

N s +2 



J I [drjo] [drj ] exp ^ // 

'Ns+2 



-Vo 



(P_ - rm q P + ) - ] [ T ; T 1 (P + - rm,P_) 

\ n n 

|[ T^P- -P + )j + l[ T s' lj - + U T s'Jr 



Vo 



>n+l 



s=l 



N s +2 



J+Jn+i J] T;\P + -rm q P_) 



s=n+l 



Ns+2 



Vo 



'Jn+l ] [ 1 P-J ~ Jn+lJn+1 



s=n+l 



Finally integrating (r)o,rjo) of (T2T1) . we obtain the generating functional 



Z\Jqi Jq-, Jni Jn+li Jn+l\ 



Ns+2 



(P_ - rm q P + ) - } [ T-\P + - rm q P-) 



Ns+2 



(24) 



(25) 



(26) 



(27) 



exp 



Jn+l ] [ T s ' P ./ + J n +lP).+ l + 



+ 



AT s +2 



j + Jn+i n 



s=n+l 



s=n+l 
Ns+2 



(P_ - rm,P + ) - JJ r7 x (P + - rm,F_ 



-i 



Ns+2 



n T *~ ip - - p +) j +n T ^ + n^ 1 ^ 

s=l s=l s=l 

Kdet[r(P c + m q )\ exp | J n+1 T v l P^ J + J n+ iJ n+1 + 

+ [J+ J n+1 T u 1 (P + -rm q P_)]r-\D c + m q )- 1 [j + f ~ x J n + f 1 J n+1 ] }, (28) 



where we have used the identity 



N s +2 



N s +2 



P+ + 1 1 T^P- P - ] [ T s '/'. 



8=1 



I±^=rP 



and defined 



iV s +2 

if 1 - n T * r 1 



s=n+l 



T- 1 = ][ PT 1 = T7 l Tv\ 



1 fr>-l 
L ' 



T- 1 = (-P+ + T _1 P_) T 
K = .7det [-P + + T- X P_] . 



Equation f[28|) is one of the main results of this paper. 

With the generating functional (|28|) . we obtain the propagators in a gauge background 
as follows. 

(I) The valence quark propagator 

1 5 2 Z 



(q(x)q(y)) 



(D c + m q ) l {x,y), 



(29) 



Z 5J q (x)5J g (y) 

where the subscript in the functional derivative denotes setting all J's to zero after differ- 
entiation. 

(II) The mixed correlator of the first kind 

(q(x)fjn(y)) = 



where 



Z 8J q (x)8J n (y) 



D^m^- 1 {-P + + T- l P_Y 1 Tl l 



D- X {m q ) = (1 + rD c )(D c + rrig)- 1 = r + (1 - rm q )(D c + m,)" 1 , 



(30) 



(31) 



the sea quark propagator. 
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(Ill) The mixed correlator of the second kind 

1 5 2 Z 



(q(x)rj n+1 (y)) 



Z SJ q (x)SJ n+1 (y) 



1 rp-1 



— (Dc + m,)- 1 {-P + + T- l P_Y l T L 

"7f D _1(m<?)75 T^+T = ^Vniy))- 



(32) 



(IV) The mixed correlator of the third kind 

1 6 2 Z 



(Vn+i(x)q(y)) 



Z 5J n+1 (x)5J q (y) 
Tv\-rm q P_ + P+)(P C + m,)" 1 ^ + T^P.y/f 



1 



T -l 



2v^ 



1 



1 + rm q 
1 — rrrir, 



75 j D~ l (m q ) 



1 — rm, 



(33) 



For completeness, we also consider the generating functional for n-point Green's function 
of fermion fields in full QCD with ODWF (satisfying the normalization condition Z[0] = 1) 

j" g— Ag — Af — Apv — ^0+j?0-P+ J— fjnJn — Jn+lVn+l— Vn+lJn+1 — VN S +3 P -J 
Z\Jqi Jq, Jfl-i Jn+l-, Jn+l] — j e -A g -A f -A PV ' 

where j = j [dU] [dip] [dip] [dtp] [d(p\ , A g is the gauge action, and Apy is the action of the 
Pauli-Villars fields {0 S ,0 S } with m q = 1/r, i.e., 

N s +3 

Apy = y^ u (f) s (x){(p s D w + T) x ,x'8s,s' + (or s D w - l) X)X i (P_6 S / )S+1 + P + 5 s / jS _i) }0 S / (x') , 

s,s'=0 x,a;' 

with boundary conditions 

P+0(x,-l) = -P + 0(x,iV s + 3), 

P_0(x, iV s + 4) = -P-</)(x, 0). 
Since the integrals over the fermion fields have been done, we proceed to evaluate the 
integrals over the Pauli-Villars fields in ( l34|) . Using the Gaussian integration formula for the 
boson fields, and following the procedures similar to above for the fermion fields, we obtain 



7T 



N s +4 Tf-1 



K- l det(l + rD c y l . 



(35) 



Substituting ([28]), and (J35J) into (J33D, we have 

1 



Z[Jqi Jq-, Jni Jn+l: Jn+l] 



j[dU]e- A 9 detP>(m q 



x 



[dU]e~ A * det D(m q ) exp \ J n +{T V L P_J + J n+1 J n+1 + 



+ [J+ Jn+iT u \P+-rm q P^)] r-\D c + m q y l [j + f' 1 J n + f' 1 J n+1 j }, (36) 
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where 

D(m q ) = {D c + m q )(l + rD c )- 1 

1/1 \ ( 37 ) 

= m q + 7} [~ ~ m i ) + 

the effective 4D lattice Dirac operator, with D c and S defined in Eqs. ©-([HI)- Setting the 
nonzero weights p s = cu s + d, and a s = cu s — d, where c and d are constants, then H s = u s H 
with H = cH w (l + djsHw) -1 . For finite N s , with the optimal weights {uj s } given in Ref. 
jsj], S is exactly equal to the Zolotarev optimal rational approximation of the sign function 
of H, i.e., S = Sopt(H) = HRz(H), where Rz(H) is the optimal rational approximation of 
(H 2 )- 1 ' 2 fl, 3- In the limit N s oo, S -»■ H(H 2 )~ 1 / 2 , and D(M is exactly equal to the 



overalap Dirac operator, satisfying the Ginsparg- Wilson relation 13] 

D(0) l5 + l5 D(0) = 2rD(0) l5 D(0). 

IV. A FORMULA FOR THE RESIDUAL MASS 

Now we are ready to derive a formula for the residual mass, in terms of the quark 
propagator in a gauge background. The denominator of (1171) can be evaluated as 

^2(q(x)X a -f 5 q(x)q(y)X b -f 5 q(y)) 

X 

= -tv F (X a X b )tv DC {[(D c + m q )- l ]\D c + m,)" 1 }^ y), (38) 

where the subscript F denotes the flavor space, and the subscript DC the Dirac and color 
spaces. In the following, the subscripts F and DC will be suppressed. 

Using Eqs. (JEZJ, ([H]), (gTJ, and (EU])-(J33]), we evaluate the numerator of (HZ]) as 

^(J 5 a (x,n)g( 2/ )A fe 75 g(y)) 
= tr(A a A b ) J ^tr[(g( 2/ )^(x))P_^ n+1 (x)g( 2 /)) 7 5] - tr[(g(y)^ + i(x))P + ^ n (x)g( 2 /)) 7 5 
= tr(A a A 6 ) <^ -J2tv[(q(y)fj n (x))j 5 P-(r] n+ i(x)q(y))j 5 ] + ^ tr[(g(y)r/ n+1 (x))7 5 P+(r/„ + i(x)g(y))7 5 

^ X X 

= tr(A a A 6 ) X)tr[( g (y)»7 n (x)>(^ + i(x)g(y)>7 5 ] 

X 

= tr(A a A 6 )^tr[ T5J D- 1 K)(r J D - P-) l5 (rD - P + )D-\m q )]{y,y) 

= -tr(A a A 6 )^tr{[ J D^ 1 K)]t(l - S 2 ) D-\m q )](y , y) . (39) 
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where D = D(0) = D c (l + rD c )~ l . Using (EH]) and ([39}, we can rewrite ([HH as 
m r „x 1 tr{[D-y g )]t(l-^)D-i(m q )}(y,y) 

mreslZ/j 4r tr{[(D c + m^]t(D c + m g )-i}(y,y) ' l4Uj 

where D~ 1 (m q ) is the sea quark propagator, and (D c +m 9 ) _1 is the valence quark propagator. 
Therefore, ( 1401) is well-defined only in the unitary limit, with the valence quark mass e qual 
to the sea quark mass. We note that Eq. ( 1401) is consistent with the form used in Refs. [U| 
and |l4j, but not in the unitary limit. 

Similarly, the global residual mass (fl8|) can be written as 

1 Tr{[D-\m q )]^l - S 2 )D-\m q )} 



M r , 



(41) 



4r 1i{[(D c + m q y i ]^D c + m q )- 1 } ' 

where Tr denotes the trace over the Dirac, color, and site indices. 

Nevertheless, it is tedious to compute the residual mass via ( 14*0]) since it involves the 
multiplication of S = (1 - fL^X 1 + IL r *) _1 = H YZ=i h i{H 2 + dj- 1 to the column 



vectors of D 



-l. 



m 



requiring conjugate gradient with multi-shift. 



In the following, we derive a practical formula for the residual mass, which only involves 
the valence quark propagator. Then the residual mass can be obtained once the valence 
quark propagator has been computed. 

We observe that the numerator of (1401) can be decomposed into two parts 



tr{[D-\m q )?D-\m q )}(y,y)-tr{[SD-\m q )}\SD-\m q ))}(y,y). 



(42) 



Using fl3~TT) and (jSJ), we obtain 



S = 75 



2r- 



D(m q ) - m q 
1 — rrrin 



- 1 



and 



SD(m q )~ l 



75 



2r 



1 + rm a 



1 — rm q 1 — rm q 



D(m n 



7 5 [r - (1 + rm q ) (D c + m q ) x ] 



Thus the second term in fT42l) can be evaluated as 



tr {[SD-\m q ))]\SD-\m q )}(y,y) 

= r 2 tr! — 2r(l + rm q )Re tr(D c + m q )~ 1 (y, y) 

+ (1 + rm q )hi{[(D c + m q )-^{D c + m,)- 1 }^, y). 



(43) 
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Using (I3TT) . the first term of (142!) is evaluated as 

ti{[D-\m q )^D-\m q )}{y,y) 

= r 2 tr! + 2r(l — rm q )Re tr(D c + m q )~ 1 (y, y) 

+ (1 - rm,) 2 tr {[(D c + m q )- l ]\D c + m^ 1 } (y, y). (44) 

Substituting (jUJ) and fj4~3l) into (HOl) . we obtain a formula for the residual mass 

Re tr{(D c + m 9 )- 1 (y,y)} 

mres ^ = + jrm Z um Z fi=T7 n ~ m <?> 45 

tr{[(D c + m 9 )T(L) c + mg)] %,?/)} 

which immediately gives the residual mass once the 12 columns of the valence quark prop- 
agator (D c + m 5 ) _1 (i,?/) have been computed. Also, it is appealing from the viewpoint 
of exact chiral symmetry, since the first term in ( I45p gives m q when D c is exactly chirally 
symmetric, (i.e. -D c 7s + 1§D C = 0), thus the residual mass is exactly zero. 
Similarly, the global residual mass (jUJ) can be written as 



Equations (145]) and P6|) are two of the main results of this paper. 

Now we consider an ensemble of gauge configurations generated in full QCD with nf 
flavors, obeying the probability distribution 

n f 

Y[detD(m f )e~ A3 , 

/=i 

then the ensemble average of the residual mass can be written as 

_ J[dU}U7 =1 detD(m f )e-^m res (y) 
" ~ RdU] TTfU det D(m,)e-+ ' 

which would be independent of y if the number of gauge configurations is sufficiently large. 
Obviously, the ensemble average of the global residual mass, (M res ), would tend to the 
limiting value with much less number of configurations. 

V. AN UPPER BOUND FOR THE RESIDUAL MASS 

For ODWF, S(H) = S opt (H), the Zolotarev optimal rational approximation of sgn(ff) = 
H(H 2 )~ l l 2 , provided that the eigenvalues of H 2 lying in the range [A 2 ^, Aj^J, where y? min 
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and y? max are the lower and upper bounds used for computing the nonzero weights {u s , s = 
1, ■ • • , n — 1, n + 2, N s + 2}. Thus, for any gauge configuration yielding eigenvalues of H 2 
lying in the range [A^ in , A^ a J, the residual mass must be bounded since it is a function of 
the sign function error 1 1 1 — S opt (H) 1 1 which is always less than dz, the maximum deviation in 
the Zolotarev optimal rational approximation. In the following, we obtain an upper bound 
for the global residual mass in lattice QCD with ODWF. 
The numerator of (HH) can be written as 



TriiDim,)- 1 }^! - S 2 )D-\m q )} = Tr{(l — S 2 )(D^D)~ 1 (m q )} 



< 



-1 



Tr{(l - S 2 )(rfD) \m q )}\ < 5>,&, (47) 



where the von Neumann's trace inequality has been used, and aj and (3j are the eigenvalues 
of |1 — S 2 \ and (Z)^) _1 (m f/ ) respectively, in the ascending order, i.e., a± < a 2 < ■ ■ ■ < dN, 
and fa < 02 < • • ■ < Pn- 

For ODWF, S = S opt , ctj < 2d z . Thus (jUD gives 

Tr{[D(m q r l )\l - S 2 opt )D-\m q )} < 2d z T,{{rf D)' 1 (m,)}, (48) 

and (jHJ) implies 

M < dz Tr{[£>(m t? )~ 1 ] t £>~ 1 (m g )} 
res ~ 2rTr{[(D c + m q )-^(D c + m q )-iy 1 1 

From (EI]), the singular values of D _1 (m (? ) and (-D c +m q ) _1 have the following relationship 

\. = r + (1 - rm q )^, (50) 

where £j is a singular value of (D c + m q )~ l and Xj is the corresponding singular value of 
Z^ _1 (m g ). Then ( 1501) gives 

2 /TWt^ -2 



£il A il 3 _,i ^< Re (0> 



r 



E,feP a ~"" s) +2r(1 - rm5) W + W)' (51) 

where 

jv JV 
^| 2 > = iV^ 12 ' ( Re (0) = ^E Re (0) ? (52) 

3=1 3=1 

and N is the total number of singular values of (D c + m q )~ l . Therefore 



TiilDim^D^img)} n , 2 . . (Re(£)) 

11 ' J "' - 1 - rm,) 2 + 2r(l - rm g ) x ;; ;/ + (53) 



j-.2 



Tr{[(D c + m q )-^(D c + m q )-i} v " " (Kl 2 > 

15 



and (jUJ) becomes 
M 



< 



dz 
2r 



;i - rm,) 2 + 2r(l - rm,) x 1 + 



<l£| 2 > ' (KI 2 )J' 

Thus, to obtain the upper bound of M res amounts to working out an upper bound for 
(Re(0)/(|f | 2 ), and a lower bound for (|£| 2 ). 

First we work out a lower bound for (|£| 2 ). The eigenvalues of V = 75^ can be expressed 
as {Rje 10 ^ ,j = 1, ■ ■ • , N}, where Rj can be bigger or less than one since the chiral symmetry 
is not exact for finite N s . Then the corresponding eigenvalues of rD c = (1 + V)(l — V)~ x 
can be expressed as 



Xj + iyj 



1-R 2 + i2Rj sin 6j 
l + R]- 2Rj cos Oj ' 



(55) 



Thus the eigenvalues of (D c + m q )~ l are rjj = r(xj + rm q + iyj)~ l . For finite N s , (\r]\ 2 ) is not 
exactly equal to (|£| 2 ), since [V"t, V] 7^ 0, due to the eigenvalues of V not falling on a circle. 
However, in estimating the lower bound of (|£| 2 ), one must fix all eigenvalues of V on a circle 
with a radius having the maximal deviation from one. Then, in this case, [V\ V] = 0, and 
2 ) = (l£| 2 )- Thus, we can use (\r]\ 2 ) to estimate the lower bound of (|£| 2 ). Using fl55l) and 



setting Rj = R, and m 



1 N 

-Y 



rm q} we obtain 

o N 



r 



3=1 
,2 



(xj + m) 2 + y 2 



-T- 

N 4-1(1 



[l + R 2 -2Rcoa9i 



1 + m 



2(1 + R 2 ) - (1 -m)(l -i? 2 ) 



+ m 2 )(l + i? 2 ) + 2m(l - i? 2 ) + 2R(1 - m 2 ) cos 9j ' 

(56) 



;i + m 2 )(l - R 2 ) + 2m(l + R 2 ) 
where we have assumed that the distribution of the eigenvalues of V is uniform in 9, and 
used the formula 

d6 1 



1 

27 



2- 



£ + A cos 9 y/B 2 - A 2 
For ODWF, |1 — R 2 \ < 2d z , and the lower bound of fl56l) is attained at R = y/1 — 2d z , i.e., 



m > 



(3 — m)d2 



(57) 



1 + m [2m + (1 — m) 2 dz 
Next we work out an upper bound for (Re(£))/(|£| 2 ). Again, using fl55l) and setting 



Ra = R, and m 



rm q , we obtain 



1 - 



r(xj + m) 

I? 



iV ' (xj + m) 2 + y 
3— 1 

1 — i? 2 



N ^ (1 



1- R 2 + m(l + R 2 - 2Rcoa9i 



;i + m 2 )(l- J R 2 ) + 2m(l + J R 2 ) 
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m 2 )(l + # 2 ) + 2m(l - i? 2 ) + 2i2(l - m 2 ) cos 0,- 
+ r- 1 m(|e| 2 ), (58 



and 



2 ) 



m + 



;i + m)(l -i? 2 



2(1 + i? 2 ) - (1 -m)(l-R 2 ) 



(59) 



where ( 156]) has been used. 



For ODWF, |1 - i? 2 | < 2d z , the upper bound of flSJD is attained at i? = VI - 2d z , i.e. 



(Mm 

(lei 2 ) 



< r _1 



m 



(1 + m)dz 
2 - (3 - m)d z 



Substituting f[6"Uj) and flBTI) into (IMl) . we obtain 



dz_ 
2r 



2(1 + m) 



2 - (3 - m)d z 



(60) 



(61) 



where m = rm q . This is one of the main results of this paper. 
In general, < dz < 0.5, this gives 

2(1 + m) 

1 + m < +^ — +^ < 4. 

V ; - 2 - (3 - m)d z ~ 

Thus the upper-bound of M res varies in the range 

— {I + m)< (M res ) upper - bound < —. 
2r r 

For ODWF, d z < 1 for A^ s > 1, then QSTJ reduces to 



C?5 



A^res < — (1 + rm a ) ~ 

- 2r v 97 2r 



(62) 



(63) 



(64) 



where m q <C mpy = r 1 has been used in the last approximation. 

Moreover, dz is an exponentially decreasing function of N s , and it can be parametrized 
6| 



as 



d z (N s , b) = A{b) exp{-C(6)iV s }, b = \ 2 max /X min , 



(65) 



where A(b) and C(b) are positive definite functions of b. This immediately implies that the 
global residual mass for lattice QCD with ODWF is an exponentially decreasing function 
of N s , regardless of the size of the lattice, at zero or finite temperatures. However, this 
scenario holds only when all eigenvalues of H 2 are falling inside the interval [A^ n , A 2 ^]. 
In general, in the simulation of full QCD with ODWF, after fixing A^ in and \ max in the 
beginning of the simulation, it could happen that some (accepted) gauge configurations in the 
course of the simulation may yield eigenvalues of H 2 lying outside the interval [X min , A^J. 
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Then the global residual mass of such "unbounded" configurations could be larger than the 
upper bound flBTj) . especially for those with many eigenvalues of H 2 less than \ 2 min . Thus, 
after generating an ensemble of gauge configurations, the ensemble averaged residual mass, 
(m res (y)) or (M res ), could be larger than the upper-bound fIBTT) . In the following, we discuss 
to what extent the low-lying eigenmodes of H 2 modify the upper-bound (I6"T|) . 

Consider a configuration U of which \H\ has N a eigenvalues (hi, % = 1, • • • ,N a ) less than 
Xmin, i.e., hi < hi < ■ ■ ■ < fiN a < Xmin- Here we assume N a <^ N = \2L Z T . For each of 
these N a eigenvalues, the corresponding eigenvalue of |1 — S 2 pt (H) \ is greater than 2d z , with 
the maximum equal to 

2d a =\l-S 2 opt (h 1 )\, (66) 
where hi is the smallest eigenvalue of \H\. Therefore (j47H is modified to 

N-N a N 

Tr{[D(m q )~ 1 ]\l - S 2 opt )D-\m q )} < 2d z £ ft + 2d a £ ft 

3=1 j=N-N a +l 

= [2(4 - d z )Q a + 2d z ]Tr{(D^Dy\m q )}, (67) 



where 



N 

Qa = ^I 1 " , (68) 



N 



^2j =N-Ng + l 



= Tz{{rfD)-\m q )}. (69) 

Now the upper- bound of Tr{(Z>t£>) 1 (m,)}/Tr{[(D c + m,)- 1 ]t(D c +m,)- 1 } can be evaluated 
as before, except replacing d z with d a . Then the upper-bound of the global residual mass 
fl6TT) is transcribed to 



d z + {da - d z )Q a 
2r 



2(1 + m) 



o -F(N.,m,N a ,h), (70) 

2 — (3 — mja 

where the factor 2(1 + m)/[2 — (3 — m)d a ] is bounded between (1 + m) and 4, similar to 
( |62l) . Thus the most significant change due to the "unbounded" low- lying eigenmodes is to 
replace d z with d z + (d a — d z )Q a , in the first factor of (I7D|) . 
Next we evaluate Q a . Using (J5D|) . we obtain 



Sg^-jy.+i ft _ (iV ffl /iV)r 2 + 2r(l - m) [Re(QJ +(l-m) 2 ne| 2 J 

Ef=i ft r2 + 2r d - ™)0MO> + (i - ™) 2 (l£l 2 > 
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where £ 3 - is a singular value of (D c + m 9 ) 1 , and 



j = N-N a + l 

N , 

(i^ 2 ) = ^Eior> (Re(0) = T7E Re fe)- 



(72) 



(73) 



j'=i i=i 
To estimate above sums, we follow the same procedure in obtaining ( 156]) and (158]) . and also 
assume that the distribution of the eigenvalues of V = 75 S opt is uniform in 9. Then we have 



(MO) = r 



2(1 + R 2 ) - (1 - m)(l - fl 2 ) 
1 + m [(1 + m 2 )(l - i? 2 ) + 2m(l + i? 2 ) 
1-R 2 



m 2 \ 

fae(0J 



1 + m 2 )(l - i? 2 ) + 2m(l + # 2 ) 

2 /-7T 



+ r- 1 m(|£| 2 >, 



1 + R 2 -2RcOs6 



7T J 9a (l + m 2 )(l + R 2 ) + 2m(l- R 2 ) + 2R(1 -m 2 )cos6'' 



d9- 



1 -R 2 + m(l + R 2 -2i?cos( 



1 + m 2 )(l + # 2 ) + 2m(l - R 2 ) + 2R(l - m 2 ) cos9'' 



(74) 
(75) 
(76) 
(77) 



where R = a/1 — 2d a , and 6 a = (1 — N a /N)7r. Then the denominator of Q a becomes 

2r 2 



r 2 + 2r(l - m)(Re(0) + (1 - ™) (10 ) 



(78) 



2m + (1 — m) 2 d a 

To evaluate the integrals (176|) and (I77j) . we perform the change of variable X — 71 ~ @-> an d 
obtain 

~ 2 l + J R 2 -2/2cos(vr-x) 



m 2 \ 

fHe(0J 



7T 



7T 





N a -K/N 



dx- 



1 + m 2 )(l + # 2 ) + 2m(l - i? 2 ) + 2i2(l - m 2 ) cos(vr - x ) 
1 - R 2 + m(l + R 2 - 2 J Rcos(vr - X )) 



,(79) 



1 + m 2 )(l + R 2 ) + 2m(l - R 2 ) + 2R(1 - m 2 ) cos(tt - X ) 



.(80) 



Since the upper limit of the integrals is N a 7c/N < 1, we can use the approximation 
cos(7r — x) — — 1 + X 2 /2 in the integrand, and obtain 

2 



nei 2 j 



rd n 



N a 



N \ m d a + (l-d a - VI -24) 



[Be(0J A " 



TV md + (1 - 4 - VI - 24) 
where we have used the formula 

A + B X 2 B {BC-AD)t^{yDjC X a) A 



(81) 



(82) 



„ dX C + D X 2 D 



Xa 



D\fCD 



c 



Xa + 0(x 2 a ) 
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and suppressed the higher order terms of 0((N a /N) 2 ). Then we obtain the numerator of 

Qa 



N, 



N 



r 2 + 2r(l - m) [Re(OJ + (1 - m) 2 [|£| 2 J 

2 



N a \ 2 f i-^T=2d a 



, r v — — . (83) 



Using QZSD and QSSD, we get 



Qa = -T7 [m + (1 - m) 



N\ 2 J \i-y/i-2d a -(l-m)d a 

N n r(l + m) + (l-m)Vl-2d a " 



L(l + m) - (1 -m)vi - 2d, 
Since < c? a < 0.5, we have 



(84) 



^O.^l- (85) 



AT - - N \m 

In other words, Q a is a monotonically decreasing function of d a , which in turn is a mono- 
tonically increasing function of N s , with the upper bound (N a /N)/m. 
For iV s > 1, d z < d a < 0.5, then 

m.-«w.=«w.=4 (^T^)- (86) 

where m < 1 has been used. If m d a /2, then (}8"6]) gives (d a — dz)Q a — 2N a /N, which is 
almost independent of N s . This immediately implies that the first factor 
[dz + (d a — dz)Q a ]/ (2r) in the upper-bound of the residual mass flTDl would look like almost 
"saturated" (decreasing slowly with respect to iV s ) after N s greater than some threshold 
value i\T* hres which depends on N a (the number of eigenvalues of \H\ smaller than A m j n ) and 
weakly on h\ (the smallest eigenvalue of \H\). In other words, if there are some eigenvalues 
of | if | smaller than A m j„, the exponential bound cannot be sustained for N s > Nl hres . This 
is one of the most interesting results emerging from our theoretical analysis. 

In Fig. [U we plot the upper-bounds (|64|) and (!70|) versus N s respectively, where the 
integrals in Q a are evaluated numerically. Here we set X min = 0.05, and X max = 6.20. The 
solid line is the upper-bound fl6"l|) of the global residual mass, for the case when all eigenvalues 
of \H\ fall inside the interval [A min , X max )- It decays exponentially with N s . The dotted 
line is the modified upper-bound (ITU]) of the global residual mass, for the case when some 
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eigenvalues of \H\ are smaller than \ m in- Here we set N a = 5 (i.e., 5 eigenvalues of \H\ smaller 
than X m in), and the smallest eigenvalue hi = 0.005 = A m j„/10. We see that the modified 
upper-bound of the global residual mass of ODWF decays exponentially up to N s ~ 18, then 
it decays like 1/A 3 for N s ~ 18-23, and almost "saturates" at ~ 10~ 5 for N s ~ 24-33. This 
agress with our theoretical analysis of the first factor [dz + (d a — d z )Q a ]/(^ r ) in the upper- 
bound of the residual mass (170|) . which would become almost "saturated" when N s > Nl hres . 
Moreover, due to the second factor 2(1 + m)/[2 — (3 — m)d a ] in (jTCTjl . the exponential bound 
dz/(2r) for N s < Nf ires is increased by a factor ~ 3, as shown in Fig. [TJ 

Next, we turn to an ensemble of gauge configurations {U{\. Let the smallest eigenvalue 
of \H\ with the gauge configuration Ui be , and the probability distribution of {h^} 
satisfies 

/ dhp(h) = 1, 

where h min = min(/?,[ 1 ' ) , hf\ ■ ■ ■ , h± \ ■ ■ ■ ) and h max = max(/i[ 1 , h\ , ■ ■ ■ , h\ , ■ ■ ■)■ Then the 
upper-bound of the global residual mass for an ensemble of gauge configurations is 

(M res ) < — U(x) f dhp(h)F(N s ,m,N a ,h) + ) [z9(z) + x0{-x)]\ ,(87) 

""max i^min \ J h min V ^ r J ) 

where x = X m in-h m in, y = min(A min , h max ), z = h max -X min , and F(N S , m, N a , h) is defined 
in (JTDjl . This is one of the main results of this paper. 

VI. NUMERICAL TESTS 

In this section, we test to what extent the upper-bound ( 18T|) of global residual mass for 
an ensemble of gauge configurations is satisfied in large-scale simulations of lattice QCD 
with ODWF. To this end, we perform hybrid Monte Carlo (HMC) simulations of two fla- 
vors QCD on the 16 3 x 32 lattice with ODWF (setting kernel H = H w ) at iV s = 16 and 
Xmin/X max = 0.05/6.20, plaquette gauge action at (3 = 6.10, and sea-quark mass m q a = 0.01. 
We have generated 2730 trajectories. After discarding the initial 300 trajectories for ther- 
malization, we sample one configuration every 10 trajectories. Thus we have 243 gauge 
configurations. For each configuration, we compute the valence quark propagator with mass 
nT-vaiU — m sea a = 0.01, and use the formula ( 145]) to obtain the residual mass. We perform the 
same calculation for 3 different cases: N s = 8, 16, 32, with the same \min/X max = 0.05/6.20. 

21 




10" 11 - 



d z /(2r) 



10- 12 



10 



15 



20 



25 



30 



N 



FIG. 1: The solid line is the theoretical upper-bound (|64p of the global residual mass, for the case 
when all eigenvalues of \H\ fall inside the interval [X m i n , X max \. Here we set \ m in = 0.05, and 
A mai = 6.20. The dotted line is the modified upper-bound (|70"j) for the case when some of the 
eigenvalues of \H\ are smaller than A m j n . Here we set iV a = 5 (i.e., 5 eigenvalues of \H\ smaller 
than X m in), an d the smallest eigenvalue h\ = 0.005. 

Then we obtain the averaged residual mass (m res ) for N s = 8, 16, 32, which are denoted by 
squares in Fig. |2j We see that for N s = 8 and N s = 16, (m res ) satisfies the exponential 
bound f )64|) . dz/(2r). However, for iV s = 32, (m res ) is much larger than the exponential 
bound. 

Next, we compute the general upper-bound ( 1H71) to see to what extent it is satisfied 
by (m res ). For each configuration, we project 250 lowest lying eigenmodes of \H W \. Then 
we obtain the smallest eigenvalue (hi) of l-H^I, and also N a , the number of eigenmodes 
of \H W \ with eigenvalue smaller than X m in- Among the set of 243 smallest eigenvalues 
{hi , i = 1, • • • , 243}, the minimum is h m i n = 6.99106 x 10~ 5 < X m i n , while the maximum is 
hmax — 0.1186 > X m i n . The probability distribution of hi can be fitted by the "log- normal" 
function 



where p = 27.4326(1.2364), a = 0.923(43), and h = 0.0108(5). The average number of 
eigenvalues of \H W \ smaller than X m i n is (N a ) = 1.778(85). The histogram of N a is plotted 
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FIG. 2: The solid line is the theoretical upper-bound dz/(2r) of the global residual mass, for the 
case when all eigenvalues of \H\ fall inside the interval [X m i n , X max ]. Here we set X m in = 0.05, 
and Xmax = 6.20. The dotted line is the general upper-bound ([87)) . The squares denote the 
average residual mass obtained with 243 gauge configurations from simulation of 2-flavors QCD 
with ODWF at (3 = 6.10, on the 16 3 x 32 lattice. 



in Fig. [3j Using ( 188)) and the information of N a , we obtain the upper-bound ([87]) of the 
global residual mass as a function of N s , which is plotted as the dotted lines in Fig. [2j We 
see that the data points of (m res ) are in good agreement with the upper-bound (JS7J). 

A salient feature emerging from this numerical study is that the exponential bound ( 164)) 
for the global residual mass holds for N s < 18. However, it is difficult to sustain the 
exponential bound for N s > 20, due to the low- lying eigenmodes of \H\ with eigenvalue less 
than \ m i n , in agreement with our theoretical analysis of the first factor [dz+(d a — dz)Q a ]/(^ r ) 
in the upper-bound of the residual mass ( 170)) . which would become almost "saturated" when 

Ng > N thres _ lg _ 20. 

At this point, it is instructive to compare the behavior of the residual mass of the ODWF 
in Fig. |5]with those of the conventional DWF. For the conventional DWF, (m res ) behaves 
ike 1/Ng for the Shamir kernel [15( , and 1 / iV s 2 for the Mobius kernel with tuned parameters 



161 ) . However, they do not possess an exponential bound for any interval of N s , unlike the 



case of ODWF. Moreover, it is straightforward to generalize the theoretical analysis in the 
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FIG. 3: The histogram of N a , the number of eigenvalues of \H W \ smaller than \ m i n , for 243 gauge 
configurations generated by HMC simulation of two flavors QCD on the 16 3 x 32 lattice with 
ODWF at N s = 16 and Xmin/^max = 0.05/6.20, plaquette gauge action at j3 = 6.10, and sea-quark 
mass m q a = 0.01. 



last section to the case of conventional DWF and show that the "saturation" phenomenon 
at large N s also holds for the conventional DWF. 



VII. CONCLUDING REMARKS 



In this paper, we have derived the axial Ward identity for lattice QCD with ODWF, by 
introducing two transparent layers at the central region of the fifth dimension, in addition 
to the two transparent layers at the boundaries for defining the quark fields [7|. From the 
axial Ward identity f in?]) , we obtain f lTB]) and f lTB 7 ]) as the local and global residual mass, for 
measuring the chiral symmetry breaking due to the finite extension in the fifth dimension. 

Since the global residual mass (|T6l) depends on the observable O, it is necessary to 
determine the residual mass of the quark for any observable O. So far, the residual mass has 
been only studied for the the pseudoscalar 0(y) = q(y)X b, j5q{y), the pion interpolator. It is 
interesting to see how the residual mass of the quark changes with respect to the physical 
observable. We have derived the generating functional for the n-point Green's function of 
fermion fields ( 1281) . which is essential for expressing the residual mass in terms of the quark 
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propagator. 

For the observable 0(y) = q{y)X b, y5q(y), we have obtained a new formula (jI5]) for the 
residual mass, which is useful in practice since it immediately gives the local residual mass 
once the 12 columns of the valence quark propagator (D c + m q )~ l (x, y) have been computed. 
For the global residual mass (Ho]) , it requires all-to-all quark propagators, or the low-lying 
eigenmodes of D = D c (l + rD c )~ l for an estimation. 

Moreover, we obtain the upper-bounds f JTOj) and (1ST)) of the global residual mass, for one 
configuration as well as an ensemble of gauge configurations in lattice QCD with ODWF. 
They provide a guideline for designing lattice QCD simulation with ODWF. That is, with 
the input values of r = [2m (l — rfmo)]" 1 and m q , how to choose the values of N s , X m i n 
and Xmax such that the residual mass meets the desired tolerance, versus the cost of the 
simulation. 

For the case when all eigenvalues of \H\ fall inside the interval [A min , A max ], only the 
second term in ( 1ST)) contributes, then M res < dz/(2r), an exponentially decreasing function 
of N s . However, if there are some eigenvalues of \H\ smaller than A m j n , then the first term 
of ( 1H71) also contributes, which makes the exponential bound only hold for N s < (N s ) thres , 
where [N s ) thres depends on X m i n and the low-lying spectrum of \H\. Moreover, the first 
term of (1ST)) also shifts the exponential bound dz/{2r) to a larger value d' z /(2r), where 
the ratio d' z /dz — 1 — 4, depending on m, \ min and the low-lying spectrum of \H\. For 
N s > (N s ) thres , the upper-bound would behave like l/N^, until it behaves like "saturated" 
(decreasing slowly with respect to N s ) at some larger N s . 

For ODWF with kernel H = H w , and plaquette gauge action with j3 = 5.95 — 6.10, then 
(N s ) thres ~ 16 - 20 for X min = 0.01 - 0.05 and \ max = 6.20. As demonstrated in Ref. Q, 
without fine tuning of X m i n , the upper bound (ISTT) gives a reliable estimate of (m res ) for 
N s = 16, with the same order of magnitude. 

The existence of a range of N s < (N s ) thres for which the exponential bound d' z / (2r) holds 
is the salient feature of ODWF, which provides a viable way to preserve the chiral symmetry 
to a good precision (say, m res a < 10~ 5 ) with a modest N s (say, N s ~ 16). 

Finally, we have a few words about the efficiency of ODWF, in comparison with other 
variants of DWF. So far, the tests in Refs. [ly] have been performed with quenched 
gauge configurations, by measuring (m res ) versus the cost of computing the valence quark 
propagators. However, in full QCD, the cost /efficiency of HMC simulation also depends on 
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the acceptance rate and the rate of topological tunnelling, which have not been addressed 



in Refs. 
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16 1 . We think it is premature to claim which lattice DWF is more efficient, 



only based on the cost of computing valence quark propagators versus the residual mass, 
without taking into account of the subtle issues (e.g., instability, topology freezing, etc.). 
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